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Abstract 

Axion fluctuations generated during inflation lead to isocurvature and non- 
Gaussian temperature fluctuations in the cosmic microwave background radi¬ 
ation. Following a previous analysis for the model independent string axion 
we consider the consequences of a measurement of these fluctuations for two 
additional string axions. We do so independent of any cosmological assump¬ 
tions. The first axion has been shown to solve the strong CP problem for most 
compactifications of the heterotic string while the second axion, which does 
not solve the strong CP problem, obeys a mass formula which is independent 
of the axion scale. We And that if gravitational waves interpreted as arising 
from inflation are observed by the PLANCK polarimetry experiment with a 
Hubble constant during inflation of Lfjnf > 10^^ GeV the existence of the first 
axion is ruled out and the second axion cannot obey the scale independent mass 
formula. In an appendix we quantitatively justify the often held assumption 
that temperature corrections to the zero temperature QCD axion mass may be 
ignored for temperatures T < Aqcd- 
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1 Introduction 


It has been known for some time now that compactihcations of string theory con¬ 
tain axion-like helds P, while it has been an open question whether these string 
axions could be the QCD axion [2]. To date, axions have never been detected, and 
it seems a practical impossibility, due to their large axion scale,^ that a string axion 
could be detected in the laboratory. However, axions give rise to isocurvature jS] 
and non-Gaussian |3] temperature fluctuations in the cosmic microwave background 
(CMB) radiation which may be measurable in the near future, even for a large axion 
scale. If future experiments are able to measure the isocurvature and non-Gaussian 
components of GMB temperature fluctuations they could bound the axion scale and 
possibily rule out string axions as the QGD axion. String axions are then one of the 
few known features of string theory that may be experimentally probed in the near 
future. 

Recently P. Fox, A. Pierce and S. Thomas (FPT) |H] analyzed the possibility that 
upcoming cosmological experiments sensitive to the polarization of the GMB (such 
as the PLANGK polarimetry experiment jH]) may be able to probe a string axion by 
measuring the isocurvature and non-Gaussian components of the GMB. For the string 
axion, FPT used the traditional model independent axion p[j, and concluded that if 
gravitational waves interpreted as arising from inflation are observed by the PLANGK 
polarimetry experiment with a Hubble constant during inflation of ifjnf > 10^^ GeV, 
then the model independent string axion cannot be the QGD axion. 

Here, we shall repeat their analysis for three different string axions. Following 
FPT we do so independent of any cosmological assumptions, allowing for the widest 
possible range of cosmological conditions within allowed constraints. Thus, our results 
are valid for any allowed cosmological scenario. For completeness, and the purpose of 
comparison, the hrst axion shall be the same axion used by FPT, the model indepen¬ 
dent (MI) axion [Tj. The second two axions, the primary concern of this paper, follow 
from the Binetruy-Gaillard-Wu (BGW) model 13IH], a class of effective supergravity 
theories derived from the weakly coupled heterotic string. These two axions are the 
BGW-QGD axion [2], which for most compactihcations of the heterotic string solves 
the strong GP problem and is therefore a legitimate candidate for the QGD axion, 
and the BGW-two condensate (BGW-2G) axion jHJ ITH] . which does not solve the 
strong GP problem, but may still have cosmological signihcance. These axions will 
be reviewed in the next section. 

Our conclusion is no different than FPT. We hnd that if gravitational waves 
interpreted as arising from inhation with a Hubble constant during inhation of iPinf > 
10^^ GeV are observed by the PLANGK polarimetry experiment, the BGW-QGD 
axion cannot be the QGD axion. This result depends critically on the fact that the 
axion scale for the BGW-QGD axion differs during inhation and condensation. We 
also hnd for a positive measurement by PLANGK that the BGW-2G axion cannot 

^One of the string axions that we will be considering, the BGW-2C axion (see section E3J, 
does not solve the strong CP problem and therefore is not a Peccei-Quinn/QCD axion. For this 
reason we will always refer to the axion coupling constant, fa, as the “axion scale” and not, say, the 
“Peccei-Quinn scale.” 
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obey a scale independent mass formula. 

In the next section we review the properties of the axions we will be considering. 
In section |S1 we present relic axion densities both with and without entropy dilution, 
the derivation of which is reviewed in Appendix El In section 0] the isocurvature and 
non-Gaussian CMB temperature fluctuations are discussed, paving the way for the 
axion constraints in section Section IHl presents our conclusions. In Appendix El 
we quantitatively justify the often held assumption that the zero temperature QCD 
axion mass may be used for temperatures T < Aqcd (an assumption that will be 
made when calculating relic axion densities in section (HI and Appendix El) • As we are 
following the analysis of FPT, the reader is referred to their excellent article [H], and 
references therein, for background and details. 

2 Axions 

We consider three different string axions, the model independent (MI) axion pQ, the 
BGW-QCD axion P and the BGW-two condensate (BGW-2G) axion jHlE]. Each 
will be described in this section, the important results being their mass and axion 
scale formulas, which are summarized in section m The MI axion was analyzed in 
jni, whose analysis we repeat in this paper for completeness and for comparison with 
the BGW axions, our primary concern. The BGW-QGD and BGW-2G axions are 
derived from the BGW supergravity effective string model mi 

We will hnd that the axion scales are very large {fa > 10^® GeV). Such large scales 
require the initial misalignment angle (see section El and Appendix El) to be extremely 
small {6i < 10“^ jH]), otherwise the relic axion density would be in contradiction with 
experiment (being larger than the observed dark matter density). Barring some 
mechanism that can set the initial misalignment angle to zero (or to an extremely 
small value), we have a highly hne tuned requirement. However, our purpose here 
does not require any such mechanism, as we are considering the axion without making 
any cosmological assumptions. We will find, in fact, that this means taking the initial 
misalignment angle to be zero, as this leads to the weakest possible bounds (see section 
El). So although there is an obvious difficulty with such large axion scales, it will not 
affect our conclusions. 

2.1 Model Independent (MI) String Axion^ 

The massless spectrum of any superstring contains the graviton, two-form gauge po¬ 
tential, dilaton and their superpartners. These helds are naturally formulated in 
terms of a linear superheld in a supergravity theory, which is dual to a formulation in 
terms of a chiral superheld ^1]. The chiral superheld contains the scalar dilaton as 
well as a pseudoscalar, the MI axion. For simplicity we present the chiral superheld 
formulation of the globally supersymmetric Lagrangian (the supergravity generaliza¬ 
tion being straightforward ^1]). The Langrangian in the chiral superheld formulation 

^In |H] the model independent string axion was referred to as the universal string axion. 
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is given by [T] 


c = -Ml j d^e ln(^ + s) + ^(^j d^esw^w^ + h.c.^ , (1) 

where S is the chiral superheld, hh“ is the superheld strength and Mp = 2.4 x 10^® 
GeV is the reduced Planck mass. Expanding the second term into component helds 
we hnd 

i J d^0SW‘W: + hx 3 -1(S + S)F;„Fr + 1(S - s)F;„Fr. (2) 

where we have used the same symbol, S, for the lowest component of the chiral 
superheld, is the Yang-Mills held strength and is the dual held 

strength. We notice the dilaton as the real part of S since it acts as the coupling 
constant for the Yang-Mills kinetic term and the axion as the imaginary part of S 
since it has the standard axion coupling to the gauge held. The QCD 6 angle is then 
given in terms of this string axion by 

e = 47r^i{S - S), (3) 


while the unihed gauge coupling constant, au, is given by 

A'jr __ 

= — = 27i{S + S), (4) 

9u 

and we take gjj 0(1). 

Expanding the kinetic term of (HD into component helds we hnd 



d^e\n{S + S) D 


Ml 2 


-d.ed>^e 

2 ^ 




where a is the axion. The axion scale, fa, is then given by 


( 5 ) 


fa/N =^ = V2^Mp ~ IQi® GeV, (6) 

where N is the axion anomaly coefficient. 

For the MI axion we use (non-stringy) mass formulas obtained from QGD. These 
are the same mass formulas used by FPT in their analysis j3]. As reviewed in Ap¬ 
pendix ^ the axion mass is temperature dependent. The zero temperature axion 
mass is calculated by considering chiral symmetry breaking in QGD (see Appendix 
ED. Onehnds 


m 


a 


2x/i fn 

-^— m. 

I + Zfa/N 


1.2 X 10"® eV 


GeV\ 

V fa/N )’ 


( 7 ) 


where z = rriu/md = 0.56, ttIt, = 135 MeV and /^r = 93 MeV is the pion decay 
constant. In Appendix El we justify the use of the zero temperature axion mass for 
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(nonzero) temperatures T < Aqcd- The finite temperature mass for T > Aqcd may 
be calculated by considering instantons at finite temperature [20] • One finds [1^] 


which is accurate for T > Aqcd, but not for too much larger.^ 


( 8 ) 


2.2 BGW-QCD String Axion 

It has been an open question whether the MI axion of the previous subsection could 
be the QCD axion. This question was recently revisited jH] in the context of the BGW 
model [HE] where it was shown that for most compactifications of the heterotic string, 
the MI axion"^ solves the strong CP problem and is therefore a legitimate candidate 
for the QCD axion. It was also found to have a modified axion scale, and being now 
model dependent, we refer to this string axion by a different name, the BCW-QCD 
axion. Our concern here is to see if cosmological observations could rule out the 
possibility of the BCW-QCD axion being the QCD axion. 

The BCW model is an effective supergravity model derived from the weakly cou¬ 
pled heterotic string, in which local supersymmetry is broken by gaugino condensation 
in a hidden sector, and is formulated in terms of the linear superfield instead of the 
chiral superfield. Whereas in the chiral superfield formulation determining which field 
is the axion is rather simple because of its coupling to the gauge fields, this is not 
the case in the linear superfield formulation, where determining the axion field can 

be subtle iHiiini. 

It is the two-form gauge potential in the linear superfield formulation that is 
related to the pseudoscalar field (the string axion) of the chiral superfield formula¬ 
tion. In the BCW model, below the condensation scale, the (unconfined) superfield 
strengths for the hidden sector gauge groups must be replaced by gaugino condensate 
chiral superfields, leading also to a replacement of the linear superfield by a vector 
superfield and a constraint condition [H- The phase of the lowest component of the 
gaugino condensate superfields is then related to the string axion. 

For a single hidden sector condensing (simple) gauge group, or for more than one 
hidden sector condensing gauge group but with identical /3-function coefficients, the 
axion is massless |H|. For the case of a single hidden sector condensing gauge group, 
however, there is an enlarged symmetry in the limit of one or more massless quarks. 
Taking into account the breaking of this symmetry was shown, in fact, to induce a 

^In the next section we will assume the axions are massless during inflation. For this to be the 
case, all contributions to the axion mass during inflation must be negligibly small. In particular, 
this includes higher dimensional operators which can contribute to the axion mass. To prevent such 
possibilities, FPT assumes that the Peccei-Quinn symmetry remains unbroken in the early universe, 
which they consider to be a mild assumption IH- For the BGW axions, the possibility of higher 
dimensional operators was considered in uniEi, where it was found, for most compactifications of 
the heterotic string, to be negligible. Thus, we do not require such an assumption for the BGW 
axions, and will therefore make no further mention of it. 

"^See footnote 121 
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mass for the axion as well as increase the axion scale compared to that of the MI 
axion [S]. For two light quarks, below the scale of supersymmetry breaking, the axion 
scale is given by 


/a/iV~ 


3 — 3 


■M„, 


( 9 ) 


where be is the /^-function coefficient for the single hidden sector condensing gauge 
group. There exists a point of enhanced symmetry at be = S/Svr^, leading to an 
increase in the axion scale and, as shown below, a decrease in the axion mass. In a 
study of electroweak symmetry breaking m the preferred range for be was found to 
be .03 < be < .04, which requires 


fa/N > 6 X 10^® GeV, 


( 10 ) 


which are the values of fa/N we will consider for the BGW-QCD axion. The preferred 
value for LSP dark matter was found to be be = .036 |T2]. The zero temperature mass 
formula was determined to be the same as dZD, 


m 


a 


1 +A/iy""' 


1.2 X 10"^^ eV 


/lO^® GeV\ 

V fa/N )■ 


( 11 ) 


As with the MI axion we use this mass formula for temperatures T < Aqcd and use 
the high temperature QGD mass formula (jHI) for T > Aqcd- 

Above the scale of supersymmetry breaking the axion scale and mass are modihed. 
If we take the supersymmetry breaking scale to be ~ TeV or above, then, as shown 
in Appendix m we will not be in need of the modihed axion mass. However, in the 
next section, we will be in need of the axion scale during inhation, which is expected 
to be well above the scale of supersymmetry breaking. In this case, we must use 


fa/N 




/2 1 
V 3 ATr%e 


Mp ~ 


10^® GeV. 


( 12 ) 


2.3 BGW-Two Condensate (BGW-2C) String Axion 

As mentioned above, for a single hidden sector condensing gauge group, or for more 
than one hidden sector condensing gauge group but with identical /3-function coeffi¬ 
cients, the axion is massless. In this case there is a nonanomalous R-symmetry. This 
symmetry is broken if there are condensing gauge groups with different /3-function 
coefficients, leading to a nonzero axion mass. Such an axion does not solve the strong 
GP problem and therefore is not a candidate for the QGD axion. It may, however, 
have cosmological signihcance, which warrants its inclusion in this paper. 

In the case of two condensates with coefficients bi > 62 , the MI axion scale is 
reduced by a factor of VQ/bi jini. As mentioned in the previous subsection, the 
preferred range of values for bi = be is .03 < bi < .04. The MI axion scale was given 
in (P) as ~ 10^® GeV, so 

/a/3V~ ^ X 10^® GeV~ 6-8 X lO^^GeV. (13) 

bi 
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For condensation scales Ai S> A 2 , the BGW-2C axion mass is approximately given 


by m 


rria 


6V2{i) 

^2 


(& 1 - 


& 2 ) 



(14) 


where i is the dilaton held in the linear multiplet formulation and ms is the gravitino 
mass which we take to be ms ~ 1 TeV. The important thing to notice about this 
formula is that it is independent of the axion scale. In the context of the weakly 
coupled heterotic string a viable scenario for local supersymmetry breaking occurs 
if the hidden sector condensation scale is given by Ai ~ 10^^ GeV HH. If we take 
{£) ~ 1 |H|, bi: ^2 ~ 10“^, 61 — 62 ~ 10“^ and A 2 > Aqcd ~ 10^ MeV we hnd 


Tiia > 8 X 10 eV. 


( 15 ) 


This contribution to the BGW-2G axion mass exists for all energies below A 2 , while 
above A 2 it disappears. 

The BGW-2G axion also receives a contribution to its mass from the QGD (zero 
temperature) axion mass formula used for both the MI and BGW-QGD axions, (|7jl 
and (HID. For the axion scale in m this gives a mass m^ ~ 6 - 8 x 10“^° eV, which 
is on the same order as (USD- Also, as shown in Appendix m the temperature at 
which this axion condenses is always on the same order as A 2 . If the condensation 
temperature were greater than A 2 then the BGW-2G axion would not obey a scale 
independent mass formula during production in the early universe. In this case, the 
cosmological consequences of this axion would not differ much from the previous two 
axions. We therefore assume that the mass formula (USD is valid for calculating relic 
densities at all temperatures and only consider the BGW-2G axion when (USD is the 
dominant contribution, and thus for masses 

rria > 10“^ eV. (16) 


We mention again that this axion does not solve the strong GP problem and is 
not a candidate for the QGD axion, but may still be cosmologically interesting, in 
particular because its mass formula (USD is scale independent. 


2.4 Summary of Formulas 

In the previous subsections we found the following masses and axion scales: 


MI; 

TJia = 1.2 X 10 ® eV 1 

GeV\ 

V fa/N )' 

fa/N ~ 10'*^ 

GeV, 

( 17 ) 

BGW-QGD: 

rua = 1.2 X 10“^^ eV 

/IQi® GeV\ 

V fa/N ) 

, /a/iV>10'® 

GeV, 

( 18 ) 

BGW-2G: 

ma ^ 10“® eV, 


fa/N ~ 10'« 

GeV, 

( 19 ) 


where the masses are to be used for temperatures T < Aqcd, except in the case of the 
BGW-2G axion where its mass is to be used for all temperatures. We see that the hrst 
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two are inversely proportional to the axion scale while the third one is independent 
of it. When deriving relic axion densities we shall therefore use the following mass 
formulas: 


m 


a 


rria 



1 GeV\ 

IJn)^ 


( 20 ) 

( 21 ) 


which we will refer to as the scale dependent and scale independent mass formulas, 
and where 


^Mi — 1.2 X lO’^ eV, ^BGW-QCD — 1.2 X lO’^ eV, Cbgw- 2 C ^10 ® eV. (22) 


For temperatures T > Aqgd we use the following mass formula for the MI and 

BGW-QCD anions: 


ma(T) ~ 2.2 X 10® eV 


/I GeV\ 

\1Jn) 


Aqgd \ / Aqcd 

200 MeV J \ T 


4 


(23) 


Above the scale of supersymmetry breaking the axion scale for the BGW-QGD 
axion is modihed. We take to be the axion scale above the scale of supersymmetry 
breaking and dehne 



(24) 


where 



jBGW-QGD ^ 


fa/N 

1018 GeV’ 


yBGW-2G ^ ^ 


(25) 


3 Relic Axion Densities 


The natural starting point when discussing axion cosmology is calculating the relic 
axion density. The details of this calculation M are by now well known (see, for 
example, iniEi), so we merely present and discuss the results here, relegating the 
details to Appendix^ (where we calculate relic axion densities for axions obeying both 
scale dependent (EDI) and scale independent m mass formulas). Special attention is 
paid to the minimal relic axion density, which will be used in section El to construct 
axion constraints. 

There are a number of different possible outcomes for the relic axion density (see 
(EZli-dZlD)- In the absence of entropy dilution, for the scale dependent axion mass 
( 1201 ) and for the axion condensing, respectively, during T < Aqcd or T > Aqcd, the 
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relic axion density is given by 


Vtah^ 


Vtah^ 


~ 1.8 X 10® 




1/2 


107 eV 

for fa/N > 2.7 X 10^^ GeV 


fa/N 
1017 GeV 

e 


3/2 




107 eV J ’ 
fa/N 


or 


" (liTiv) (ira) + 

for fJN < 2 X 10'^ GeV. 


( 26 ) 


( 27 ) 


while for the scale independent mass m, still in the absence of entropy dilntion and 
for any condensation temperatnre, 


Qr.h'^ ~ 3.5 X 10^ 


C /‘-'V fJN \ 


10-9 eV 




In the presence of entropy dilntion, for the scale dependent (PUj) or scale independent 
m axion mass, the relic axion density is given by, respectively. 




for fa/N < 3.0 X 10^9 GeV 

^RH 


or C > 3.3 X 10"^^ eV 


107 eV 
2 


6 MeVy 


T, 


RH 


6 MeV 


( 29 ) 

( 30 ) 


In f|20|l - pn|l h = ifo/(100 km s~^ Mpc“^), V is a model dependent nnmber on the 
order of one, f{0f) is a correction for anharmonic effects of the axion potential (which 
for small 6i is ronghly one), Trh is the reheating temperatnre for the late decaying 
particle that dilntes the axion condensate (and must satisfy Trh ^ 6 MeV), (6**)^ + 
(Tg = (Of) is the mean square of the initial misalignment angle, 6i, averaged over the 
universe, and we have used Aqcd = 200 MeV. 

The relic axion densities are all written in terms of the quantity 

i^i) = ( 31 ) 


In the absence of inflation, each causally connected region of the universe would have 
a different (randomly selected) initial misalignment angle in the range —tt < < tt, 
and the rms value would be (Of) = 7r^/3. In an inflationary universe, however, 
string axions exist before inflation and so the observable universe would have a single, 
constant, initial misalignment angle, since there existed a time in the early universe 
when the entire observable universe was in causal contact. In this case, to use the 
rms value would be to average over each inflationary region, which tells us nothing 
about our observable universe. Thus, it would not be correct to use the rms value for 
an inflationary universe. 
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Since we are assuming an inflationary universe we now interpret all quantities to 
refer to the observable universe only. Thus, 6i, the initial misalignment angle, takes 
a constant value over the entire observable universe, so {9i) = 6i. ae is the standard 
deviation of the initial misalignment angle, but since the observable universe has a 
single initial misalignment angle, it should be zero. Below we shall reinterpret ag as 
a characterization of quantum fluctuations, which are nonzero, and so shall leave it 
in equations. 

Axions induce isocurvature [3] and non-Gaussian j3] components in CMB tem¬ 
perature fluctuations. These effects, as explained in section 01 may be parametrized 
in terms of the square and cube of the relic axion density, and the weakest bounds 
will be given by the minimal relic axion density, which itself is given by zero initial 
misalignment angle. For zero initial misalignment angle, the relic axion densities are 
determined purely by quantum fluctuations during inflation (characterized by ag). 

All massless helds undergo such quantum (de Sitter) fluctuations during inflation, 
and as long as the axion is (effectively) massless in the early universe,^ it will too 
jl8j . The quantum fluctuations can be written in terms of the Hubble constant during 
inflation, ifinf, as [TH] 

- - ^ - /‘TO- 

where is the axion scale (above the scale of supersymmetry breaking and therefore) 
during inflation and fg = fa/fa- The current bound on is TTinf < 2.8 x 10^^ GeV 

H- 

As mentioned above, we will be in need of the minimal relic axion densities. 
From (Eni)-® we see that this corresponds to vanishing initial misalignment angle. 
Replacing ag in favor of Ffinf we hnd for the minimal relic axion densities, in the 
absence of entropy dilution, for the scale dependent axion mass (HD and for the 
axion condensing, respectively, during T < Aqcd or T > Aqcd, 


~ 4.6 X 10-® /2 

Cl ^ i 


~ 6.2 X 10"^ /, 


10^ eV 

for fa/N > 2.7 X 10^^ GeV 

e 


^ y/2 ^1017 GeV 






-f 

1013 GeV ) 


-5 f2 

" ' W eV 

15 


W eV ' ’ 

GeVy/® 

fa/N ) GeV 


H. 


inf 


for fa/N < 2 X GeV, 


(33) 


(34) 


while for the scale independent mass (ED), still in the absence of entropy dilution and 
for any condensation temperature. 


~ 8.8 X 10 


-5 


I 


c 


10-9 eV 


1/2 


/ i^inf y 

V1013 Gevy 


®See footnote 01 


(35) 
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In the presence of entropy dilution, in which the minimal relic axion density is given 
for maximum entropy dilution, ^RH — 6 MeV, for the scale dependent (unD or scale 
independent eu axion mass, the relic axion density is given by, respectively. 


^mm^2 _ g_4 ^ 10 - 

d ^ s 


( i/inf y 

V1013 Gevy 


for fa/N < 3.0 X 10^0 GeV 
OT C > 3.3 X 10“^'^ eV. 


These will be used in section |S1 to construct axion constraints. 


(36) 

(37) 


4 CMB Temperature Fluctuations 


Quantum, spatially-dependent perturbations of helds during inflation can lead to 
observable temperature fluctuations in the CMB. For example, inflationary pertur¬ 
bations of the inflaton, which is an example of an adiabatic perturbation, leads to 
adiabatic fluctuations in the CMB. Adiabatic perturbations are characterized by fluc¬ 
tuations in the total energy density, Sp ^ 0, without fluctuations in the local equation 
of state, 5{ni/s) = 0. Our concern, however, will be with inflationary perturbations 
of the axion, which is an example of an isocurvature perturbation, and leads to 
isocurvature fluctuations in the CMB jH]. Isocurvature perturbations, characterized 
oppositely to adiabatic perturbations, correspond to fluctuations in the equation of 
state, 5{ni/s) = 0, without fluctuations in the energy density, 5p ^ 0. It is these 
isocurvature fluctuations, which contain the non-Gaussian fluctuations, that will lead 
to axion constraints. 

We assume that the axion is the only contributor to isocurvature perturbations 
of the CMB (all other helds undergoing adiabatic perturbations) since this gives the 
weakest bounds. Assuming otherwise would strengthen bounds by decreasing the 
axion isocurvature contribution for a (hxed) measured total isocurvature component 
of the CMB. 

We repeat here the derivation found in jS] (to which the reader may look for 
details). Isocurvature huctuations of the axion may be characterized by 


^ 5{na/s) 

Ua/s 


-3 — , s oc T , 

rin, I 


(38) 


where we shall dehne exactly what we mean by 6 below. By considering the huctu- 
ation of the total energy density, it may be shown that initially, just after the axion 
condensate is formed, the isocurvature temperature huctuations, ST/T in (IH^ . satisfy 
(5T/T)init -C {Sna/na)init- As na/s, and hence Sa, is conserved, it follows that 





(39) 
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Sa may now be written in terms of the axion field, or equivalently, the 6 an¬ 
gle. Since we are assuming an inflationary universe, we consider only our observable 
universe. In this case, the initial misalignment angle, 6i, is a constant (over the 
observable universe) and we dehne 69 to be the spatially-dependent quantum fluctu¬ 
ation and 6 = 6i + 66. Brackets () will denote averaging over the observable universe, 
so {9i) = 9i since it is constant, and we take the fluctuation, 66, to be Gaussian 
distributed with zero mean, so {66) = 0 . 

We are now in a position to dehne what we mean by 6 in (Id 8 |l and (P|). For 
concreteness, we do so specihcally for 6na. Since we are dealing with spatially depen¬ 
dent quantum huctuations of the axion, and thus of the axion number density, Ua is 
taken to be the average value of the axion number density (and is the actual number 
density in the limit of vanishing quantum huctuations), averaged over the observable 
universe. 6na is then dehned as the diherence between the actual number density at 
some point in space and the average value, Ua- To express these dehntions in terms 
of the 6 angle, we note from (EDI) that Ha oc ( 6 *^), and so 


6na oc 6*^ — (0^) . (40) 

Along with ( 02 ) = ^ 02 ) ^ _ q 2 ^ ^2 _ ^^ 2 ) _ ^ 0)2 _ ^^^^) 2 ) 

as a characterization of the average quantum huctuation of 0 , (innD may be written in 
terms of the 0 angle as 


Sa^ 


6nn 


Ua 


{6i + 66f - ( 02 ) _ 26,66 + (50)2 - al 


init 


( 02 ) 


from which we shall need 

^ 2a2 


2 20.^ + 


+ ^e) 


2',2’ 




61 + a2 


.4 30 ? + al 




(41) 


(42) 


where we have used {(69)^) = 3a^ and ((50)®) = 15a^, following from 50 being 
Gaussian distributed. 

Finally, the GMB temperature huctuations induced by isocurvature perturbations 
from an axion are given by (see jS] for details) 


f^) 6 6 G, 26,66 + (50)2 - 


where is the density of all non-relativistic matter and Qa is the relic axion density 
presented in the previous section. This equation takes into account Sachs-Wolfe con¬ 
tributions IZQ which are red shifting perturbations by matter to the GMB photons 
while they are on their way to the detector. From WMAP = TSSI qqq [213 

a AGDM cosmology. 


Isocurvature 

The observable ehects of isocurvature huctuations are conveniently parametrized 

^ ^ (mT)L} ^ , fey , 202 + a2 

((5T/T)L) Visj ((5T/T)L) '^(0? + a2)2’ 
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where COBE has measured {{ST/T)‘l^^) ' ~ 1.1 x 10“^ (22] and a current conservative 


bound for a is a < .4 1221 - 

In the next section we will rewrite this equation in terms of iJinf and the mini¬ 
mal relic axion densities of the previous section, thereby obtaining a lower bound on a. 


Non-Gaussianity 


Although the axion experiences Gaussian fluctuations {SO was taken to be Gaus¬ 
sian), it induces both Gaussian and non-Gaussian fluctuations in the GMB j3]. The 
hrst term on the right hand side of il, proportional to 26iS6, is obviously Gaussian 
distributed, but the remaining term, proportional to {S6)^ — <7 is not. This is the 
axion induced non-Gaussian fluctuation in the GMB. 

Gaussian fluctuations with zero mean are completely defined in terms of two point 
functions since all odd point functions vanish and all even point functions are products 
of two point functions. Thus, a non-vanishing three point function may act as a 
measurement of non-Gaussianity. The observable effects of non-Gaussian fluctuations 
induced by axions are then conveniently parametrized by the dimensionless skewness. 



Note that we have assumed that the axion is the only source of non-Gaussianity,® 
as this leads to the weakest bounds. Assuming otherwise would strengthen bounds 
by decreasing the axion contribution to non-Gaussianity for a (hxed) measured total 
non-Gaussian component of the GMB. GOBE has measured —S's iso < -06 [22] • 

In the next section we will rewrite this equation in terms of iJjnf and the minimal 
relic axion densities of the previous section, thereby obtaining a lower bound on 



5 Axion Constraints 

We may now use our results to constrain the axion. As mentioned in the introduction, 
we shall do so under the minimal possible assumptions, leading to the weakest possible 
bounds. 

The axion induces isocurvature and non-Gaussian components in GMB tempera¬ 
ture fluctuations as described by dH and (021), respectively. A look at these equa¬ 
tions tells us that a and —S^ are minimized for minimal relic axion density, G™™. We 
found in section |21 that minimal relic axion density corresponds to vanishing initial 
misalignment angle. Thus we may use the fact that 



( 46 ) 


®In certain scenarios, the axion cannot be the dominant source of non-Gaussianity |22| . This 
does not affect our conclusions since additional sources could only strengthen bounds. 
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to constrain the axion. The minimal relic axion densities were given in (ISl-dSZl). 
Plugging these in and rearranging, we hnd, in the absence of entropy dilution, for the 
scale dependent axion mass (isni) and for the axion condensing during T < Aqcd, 


fa/N > 3-1 X 10 ^^ 


falN> 


a 

6.1 X 10^9 GeV 




-^3,iso)2/3 


ft 


107 eV 


H, 


inf 


1013 GeV 


4^ 


H, 


inf 


107 eV / V 1013 GeV 


while for the axion condensing during T > Aqcd, 


for fa/N > 

2.7 X 10^^ GeV 




107 eV 
(47) 


falN> 


fa/N> 


4.4x1018 GeV 12 / 5 / ^ \ 

a3/5 Vl07eV; VlO^^GeV; 

6.6 X 1018 GeV , 2/5 ( ^ ( ^inf \ 

(-%so)2/5 \WeVj VlO'^GeVy* 


12/5 

12/5 


for fa/N < 

2 X 10^® GeV. 


(48) 

Still in the absence of entropy dilution, but for the scale independent mass m and 
for the axion condensing at any temperature. 


3 4^4 


a > 1.1 X 103/^ 

-^3,iso > 1.2 X 


c 




inf 


10-9 eV / V 1013 GeV / ’ 


VlO'^GeVt 


( 49 ) 


In the presence of entropy dilution, axions obeying respectively the scale dependent 
dm) or scale independent m axion mass formulas have identical constraints, but 
with different domains of validity. 


a > .10 f 
~ 1013 GeV 


-53,iso>.093/.M^^ 


for fa/N < 3.0 X 10^° GeV 




107 eV 


or y > 3.3 X 10 eV. 


(50) 

(51) 


As mentioned in Appendix 0 in the case of entropy dilution, but outside the domains 
dm or m, the above formulas are no longer valid, and the axion is constrained by 

fiTfji or dm- 


BGW-QCD String Axion 

The constraints listed above are for generic axions satisfying the scale dependent 
dm and scale independent CD axion mass formulas. Each constraint contains factors 
of fs, where 7 ^ 1 means that the axion scales during inflation and condensation 
differ. For the three axions we are considering, only the BGW-QGD axion has fs 7 ^ 1, 
as seen from dm- We now construct the specihc BGW-QGD constraints by plugging 
in its values of y and fs- 
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In the absence of entropy dilution and for the BGW-QCD axion condensing during 
T < Aqcd, 


fa/N < (3.0 X 10'^ GeV) 


4/3 


for fa/N > 3.2 X 10^^ GeV 


fa/N < (2.4 X 10'^ GeV) 


1013 GeV\ 




inf 


■) J 


(52) 


while for the BGW-QCD axion condensing during T > Aqcd, 

12/7 


/I 013 HpVX 

fa/N < (3.0 X lOi^ GeV) [ -] 

V Ainf / 


for fa/N < 2 X 1013 GeV, 


12/7 


fa/N < (2.2 X lOi^ GeV) (-^3 




V H, 


(53) 


while in the presence of entropy dilution, 

'1013 GeV 


/ ^ 

fa/N < (1.8 X I 013 GeV) Cfi/^ i - 


H: 


inf 


fa/N < (1.5 X 1013 GeV) 


1013 GeV 


I 


for fa/N < 3.6 X 10"ii GeV 


ifinf / 


(54) 

As mentioned above and in Appendix in the case of entropy dilution, but outside 
its domain, dSl is no longer valid, and the axion is constrained by (ESI- 


6 Conclusions 

The axion constraints of the previous section depend on the Hubble constant during 
inflation, hfinf. The PLANCK polarimetry experiment is expected to be sensitive to 
values of iLinf much smaller than the current bound of ifinf < 2.8 x IQi^ GeV ng. 
After briefly describing the PLANCK polarimetry experiment and its expected level 
of sensitivity, the constraints of the previous section will be examined under of the 
possibility of a precise determination of by PLANCK. 

As mentioned in section El all massless helds undergo quantum fluctuations during 
inflation, which includes the metric. Metric fluctuations generated during inflation 
appear as gravitational waves which induce tensor B-modes in CMB temperature 
fluctuations which PLANCK will attempt to measure. Just as axion fluctuations 
can be written in terms of iPmf, so too can metric fluctuations. Thus, a positive 
measurement by PLANCK of gravitational waves interpreted as arising from inflation 
would be a measurement of iPinf. A reasonable estimate is that PLANCK will probe 
H\ni > 10^3 GeV PHE]. In hgures UHSl the constraints of the previous section are 
plotted. In these hgures, the PLANCK polarimetry experiment is assumed to probe 
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the parameter space to the right of the dashed vertical lines, which are given by 

//inf ^ 1013 GeV 


ai/4 ~ (.4)1/4 

//inf . 1013 GeV 


> 


(-53,iso)l/l* ~ (.06)1/6 


= 1.3 X 10i3 GeV, 
= 1.6 X 10i3 GeV. 


(55) 

(56) 


Figure HI presents the constraints (jUj), and (IKnjl for an axion obeying the 
scale dependent QGD mass formula 0 , and with /,, = 1 , which means that the axion 
scale during inflation and condensation is the same. This includes the MI axion (but 
not the BGW-QGD axion which has fs 7 ^ 1) and these plots are equivalent to the 
plots in |H] . In the absence of entropy dilution, the allowed parameter space is above 
the diagonal lines (and we may ignore all other lines in the hgure). The diagonal 
lines do not connect because they are valid in disconnected regions. The constraints 
between them should presumably be given by a smooth interpolation which connects 
them. Such an interpolation is given by the dots. The allowed parameter space, in 
the absence of entropy dilution, is then above the dots. In the presence of entropy 
dilution we ignore the diagonal lines and dots, and the allowed parameter space is to 
the left of the thick vertical line. The thick vertical line stops because the constraint 
it describes is only valid in the region given by CT. Outside this region, and thus 
above where the thick vertical line stops, it is the diagonal lines (and dots) which are 
the proper constraint, even in the presence of entropy dilution. 

In the absence of entropy dilution, we can see in hgure Q that if PLANGK makes 
a positive measurement of gravitational waves interpreted as arising from inhation, 
and thus we are somewhere to the right of the dashed vertical line, then fa/N > 10 ^° 
GeV (so as to be above the dots). This rules out the existence of the MI axion, 
which requires, from O. /a/iV 10^6 GeV. In the presence of entropy dilution, 
hgure ^a) shows a sliver of parameter space (the space between the vertical lines) in 
which a positive measurement by PLANGK could be consistent with the MI axion. 
However, this possibility is ruled out in hgure ^b). Thus we hnd that for a positive 
measurement of gravitational waves interpreted as arising from inhation by PLANGK, 
the existence of the MI axion is ruled out jH]- 

Figure 121 presents the constraints for the BGW-QGD axion. The BGW- 

QGD axion obeys the scale dependent mass formula (ED and has fs given by (USD) 
which means that the axion scale during inhation dihers from that during conden¬ 
sation. In the absence of entropy dilution, the allowed parameter space is below the 
thin diagonal lines (and we may ignore all other lines in the hgure). The thin diagonal 
lines do not connect because they are valid in disconnected regions. The constraints 
between them should presumably be given by a smooth interpolation which connects 
them. Such an interpolation is given by the dots. The allowed parameter space, in 
the absence of entropy dilution, is then below the dots. In the presence of entropy 
dilution we ignore the thin diagonal lines and dots, and the allowed parameter space 
is below the thick diagonal line. The thick diagonal line stops because the constraint 
it describes is only valid in the region given by ( 1 ^ . Outside this region, and thus 
beyond where the thick diagonal line stops, it is the thin diagonal lines (and dots) 
which are the proper constaint, even in the presence of entropy dilution. 
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(GeV) 



(GeV) 


Figure 1: Constraints for axions obeying the scale dependent QCD mass formula 
((EDI) with the values in (j22j) ) and with the same axion scale during inflation and 
condensation. This includes the model independent string axion and these plots are 
equivalent to the plots in jH] . (a) is for the isocurvature constraints while (b) is for the 
non-Gaussian constraints. The diagonal lines correspond to constraints in the absence 
of entropy dilution, (gZl) and (HHl) . The dots are a smooth interpolation connecting 
them. The allowed parameter space is above. The thick vertical line corresponds to 
constraints in the presence of entropy dilution (H), the allowed region being to the 
left. The space to the right of the dashed vertical line is a reasonable estimate of the 
parameter space that the PLANCK polarimetry experiment will probe. 
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(GeV) 


Figure 2: Constraints specifically for the BGW-QCD axion. (a) is for the isocurvature 
constraints while (b) is for the non-Gaussian constraints. The thin diagonal lines 
correspond to constraints in the absence of entropy dilution, (El and (1^ . The dots 
are a smooth interpolation connecting them. The allowed parameter space is below. 
The thick diagonal line corresponds to constraints in the presence of entropy dilution 
(El, the allowed region being below. The space to the right of the dashed vertical 
line is a reasonable estimate of the parameter space that the PLANCK polarimetry 
experiment will probe. 
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In the absence of entropy dilntion, we can see in figure |21 that if PLANCK makes 
a positive measurement of gravitational waves interpreted as arising from inflation, 
and thus we are somewhere to the right of the dashed vertical line, then fa/N < 10 ^”^ 
GeV (so as to be below the dots). This rules out the existence of the BGW-QCD 
axion, which requires, from (uni), fa/N > 6 X 10^^ GeV. In the presence of entropy 
dilution, if PLANGK makes a positive measurement then fa/N < 10^® GeV (so as to 
be below the dots). Again, the existence of the BGW-QGD axion is ruled out. Thus 
we hnd that for a positive measurement of gravitational waves interpreted as arising 
from inflation by PLANGK, the existence of the BGW-QGD axion is ruled out. 

Figure El presents the constraints (jm|) and m for an axion obeying the scale 
independent mass formula m and with /^ = 1, which means that the axion scale 
during inflation and condensation is the same. This includes the BGW-2G axion. In 
the absence of entropy dilution the allowed parameter space is below the diagonal line 
(and we may ignore all other lines in the figure). In the presence of entropy dilution 
we ignore the diagonal line and the allowed parameter space is to the left of the thick 
vertical line. The thick vertical line stops because the constraint it describes is only 
valid in the region given by m- Outside this region, and thus below where the thick 
vertical line stops, it is the diagonal line which is the proper constraint, even in the 
presence of entropy dilution. 

In the absence of entropy dilution, we can see in hgureElthat if PLANGK makes 
a positive measurement of gravitational waves interpreted as arising from inflation, 
and thus we are somewhere to the right of the dashed vertical line, then ( < 10“^^ eV 
(so as to be below the diagonal line). This rules out the possibility of the BGW-2G 
axion obeying a scale independent mass formula, which requires, from (dni), C = ma> 
10“® eV. In the presence of entropy dilution, hgureEKa) shows a sliver of parameter 
space (the space between the vertical lines) in which a positive measurement by 
PLANGK could be consistent with the BGW-2G axion obeying a scale independent 
mass formula. However, this possibility is ruled out in hgureEKb). Thus we hnd that 
for a positive measurement of gravitational waves interpreted as arising from inhation 
by PLANGK, the possibility of the BGW-2G axion obeying a scale independent mass 
formula is ruled out. 

Just as in we have shown that if the PLANGK polarimetry experiment makes a 
positive measurement of gravitational waves interpreted as arising from inhation with 
a Hubble constant during inhation of iLinf > 10^^ GeV, then a QGD axion satisfying 
the usual mass formulas, (|7j) and (jHI), must have fa/N > 10^° GeV. In QGD, any axion 
scale and corresponding axion mass solves the strong GP problem, so this result has 
no direct bearing on the status of the QGD axion. However, it is interesting to ask 
if there are theories which can supply such large axion scales. It is well known that 
string axions have naturally large scales, so we have asked in this paper, following [^, 
whether two diherent string axions could be the QGD axion. We have shown that if 
the PLANK polarimetry experiment makes a positive measurement of gravitational 
waves interpreted as arising from inhation, then neither the model independent axion 

nor the BGW-QGD axion^ can be the QGD axion, their existence being ruled out. 

^We note that for the BGW-QCD axion this result depends critically on the fact that the axion 
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ff,nf/a'^^(GeV) 



(GeV) 


Figure 3: Constraints for axions obeying the scale independent mass formula m- 
This includes the BGW-two condensate axion. (a) is for the isocurvature constraints 
while (b) is for the non-Gaussian constraints. The diagonal line corresponds to con¬ 
straints in the absence of entropy dilution (m, the allowed region being below. The 
solid vertical line corresponds to constraints in the presence of entropy dilution m, 
the allowed region being to the left. The space to the right of the dashed vertical 
line is a reasonable estimate of the parameter space that the PLANCK polarimetry 
experiment will probe. 
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We have also shown that a third string axion, the BGW-two condensate axion, which 
does not solve the strong CP problem and is therefore not a candidate for the QCD 
axion, cannot obey a scale independent mass formnla, snch as (fT6|l . if PLANCK makes 
a positive measnrement of gravitational waves interpreted as arising from inflation. 
Finally, in Appendix El we have given a qnantitative calculation of the corrections to 
the zero temperature QCD axion mass in the context of the chiral Lagrangian and 
justihed the often held assumption that such corrections are negligible. 
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A Axion Production 

In this appendix we review the production of axions in the early universe and the 
resulting relic axion densities. We do so for the standard axion mass formulas as well 
as a scale independent formula. 

There are three methods in which axions are produced in the early universe (see, 
for example, dl): (1) thermal interactions, (2) decay of axionic strings, and (3) 
relaxation of an initially misaligned 6 angle. For axion scales, fa, much above 10^° 
CeV, which we shall certainly be assuming, axions interact so weakly that a thermal 
population never results. In an inflationary universe, which we shall also be assuming, 
strings axions exist before and during inflation, and so all axionic strings are inflated 
away.® This leaves only the third method uni for axion production, which we now 
review. 

The equation of motion for a homogeneous axion, a, which is related to the 6 angle 
by a = {fa/N)6, where N is the axion anomaly coefficient, in an FRW cosmology is 

d + 3Ha + m‘l{T)a = 0, (57) 

where H is the Hubble parameter and we have assumed that a is small in writing the 
last term (which will be corrected for below). It is important to note that the axion 
mass, rriaiT), is temperature dependent. 

At high temperatures the axion mass is very small (the axion mass formulas were 
presented in section E and will be written again below) and the mass term in (157j) 
can be ignored. If we assume initially that a = 0 then the equations of motion 

scale during inflation is given by m and not ©• 

®I am thankful to A. Vilenkin for bringing a misunderstanding of mine on this issue to my 
attention. 
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are solved by a constant axion field, a* = {fa/N)6i, 6i being known as the initial 
misalignment angle. As the temperatnre falls the axion mass increases. Once we 
reach the oscillation temperatnre (also known as the condensation temperatnre), Tqsc, 
defined by 3H{Tosc) = maiTosc), the axion begins to oscillate aronnd its vacnnm valne, 
and is said to have condensed. 

From a look at (EZD, the energy density is given by 


Pa = + ^ml{T)a‘^. (58) 

It is not hard to show that even while the axion mass is decreasing, the nnmber density, 
Ua, scales like non-relativistic matter (see, for example, El). Assnming there is no 
entropy dilntion, this allows ns to calcnlate the cnrrent axion nnmber density. 


PaO 

^aO 


iria S 


(59) 


where ma = rriaiT = 0) is the zero temperatnre axion mass and we evalnate Ua/s at 
the moment the axion condensate forms, so 

(9?) /(9?), (60) 

which follows from (EHD, along with fc — 1.44, which is a nnmerical correction for 
the temperatnre dependence of the axion mass and f{0f) which is a correction for 
anharmonic effects of the axion potential (and corrects for onr assnmption of small a 
in writing the last term in (EZI), which for small 6 is ronghly eqnal to one Hi- 
is the mean sqnare valne of the initial misalignment angle, discnssed in section El 
The final piece of information that we need is the axion mass. Axion masses for 
the three string axions we are considering were presented in section |21 For two of the 
axions (the MI and BGW-QCD axions) the zero temperatnre masses were inversely 
proportional to the axion scale, while for the third axion (the BGW-2C axion) the 
mass was independent of the axion scale. Gonseqnently, we will derive relic axion 
densities for the following two mass formnlas. 


m 


a 


m 


a 



1 GeV'i 

TWJ’ 


(61) 

(62) 


which we will refer to as the scale dependent and scale independent mass formnlas, and 
where ^ and ( are constants (which have dimensions of mass). The scale dependent 
mass (inii) will be nsed for temperatnres T < Aqcd, which we jnstify in Appendix 
El while the scale independent mass (P|) will be used for all temperatures, which we 
justify below. For axions which obey the scale dependent formula (iniD, but condense 
at temperatures T > Aqcd, we will need the high temperature axion mass formula, 
which may be calculated by considering instantons at finite temperature m- One 
finds [IS] 


ma{T) ~ 2.2 X 10^ eV 


/I GeV\ 


f Aqcd \ / Aqcd 

V200 MeVy V T 


4 


(63) 


21 






which is accurate for T > Aqcd, but not for too much larger. 

With these mass formulas, 3H = ma{Tosc), the Friedman equation, p = 
and the thermal energy density in a radiation dominated era, p = 
where = 61.75 for Tqsc just above Aqcd and g^ = 10.75 for Tosc just below Aqcd, 
we may solve for the oscillation temperature. For axions obeying the scale dependent 
mass formula (|F)T|l which condense, respectively, during Tqsc < Aqcd or Tqsc > Aqcd, 


Tosc ^ 3.3 X 10^ MeV 
Tosc — 7.0 X 10"^ MeV 


e y/viGevy/^ 

T^J \JJn) ’ 

^ Aqcd ( I GeV y/^ 

,200 MeVy* V fa/N ) ’ 


(64) 

(65) 


while for axions which obey the scale independent mass formula (P|) and condense 
at any temperature,® 


Tosc — 2 X lO’^ MeV 



( 66 ) 


We may now form the relic axion density, Gq = Pao/Pc where Pc = 2hlMp is 
the critical density and Mp = 2.4 x 10^® GeV is the reduced Planck mass. Using 
s = (27rV45)p*s7;tc) 9*s = 61.75 for Tosc just above Aqcd and p„ = 10.75 for Tosc just 
below Aqcd, with Tosc given above, and Sq = (27rV45)p*so^o^, g*so = 3.91, Tq = 2.73 
K, we have for axions obeying the scale dependent mass formula dnn) which condense, 
respectively, during Tosc < Aqcd or Tosc > Aqcd, 

-1.8 X10-" (^) m mr, (67) 

^ 3.6 X 10-^^ (r^) ' <"?) (68) 

while for axions which obey the scale independent mass formula (P|) and condense 
at any temperature, 

cx 1.1 X 10-^‘ («?> /(ft)t (69) 

where h = Ho/{100 km s“^ Mpc“^) and we have taken Aqcd = 200 MeV. Using (IMI) 
and (insi) the domains of validity for (inzD and (El may be rewritten as 

Tosc < Aqcd ^ fjN > 2.7 x 10^° GeV , (70) 

Tosc > Aqcd ^ /a/V < 2 x 10^^ GeV, (71) 

depends on the value of Tosc and in this sense this formula is not valid for all temperatures. 
However, T^sc 9* and the error in using 5 * = 75.75 (which we use in this and future equations 
for the scale independent mass El) for all oscillation temperatures is small. 
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where again we have taken Aqcd = 200 MeV. 

In section 12.41 the numerical values for the axion masses and scales we will be using 
in this paper were given. With these values, dzni) and CB tell us that for axions 
obeying the scale dependent mass formula Tqsc ~ Aqcd and below, and that 
in some cases we are in fact considering oscillation temperatures where neither the 
high temperature nor the zero temperature mass formulas are accurate. Presumably, 
in these cases, the correct relic axion density is given by the smooth interpolation 
between the densities derived from high and low temperature mass formulas. We also 
see that we will never be considering oscillation temperatures above a TeV. One of 
the axions we are considering (the BGW-QCD axion) has its mass and scale modified 
above the scale of supersymmetry breaking (see section . If we take this scale to 
be ~ TeV or above then we will not have to take this into account here (we will have 
to take this in account in section El when we consider axion fluctuations generated 
during inflation). For the axion obeying the scale independent mass formula (j(i2|l . we 
hnd that Tosc ~ A 2 , where above A 2 the scale independent mass formula is no longer 
valid. As mentioned in section (231 we will assume the scale independent mass formula 
is valid for all oscillation temperatures since otherwise the cosmological consequences 
of this axion would not differ much from the others. 

So far everything has been derived in the absence of entropy dilution. We now 
consider the possibly of a late entropy release diluting the axion condensate by 
the late decay of a massive, non-relativistic particle. For the late decaying particle to 
dilute the axion condensate it must decay after the axion condensate is formed, and 
therefore have a reheating temperature satisfying Trh ^ ^osc- However, any entropy 
release is constrained by the success of big bang nucleosynthesis, which in this case 
requires Trh ^ 6 MeV To calculate the relic axion density we must form Ua/s 
after the particle has decayed. The result is 



( 72 ) 


Using (jM|l and we may rewrite the domain of (f7^ . Tqsc > Trh, for axions obeying, 
respectively, the scale dependent dnn) and scale independent (EH) mass formulas, as 



( 73 ) 


( 74 ) 


Outside these domains, the axion condenses after the particle decays, which is during 
a radiation dominated era, and therefore the axion density is given by (EH) or El- 

B Low Temperature QCD Axion Mass 

In this appendix we explicitly calculate temperature corrections to the zero temper¬ 
ature QCD axion mass using the chiral Lagrangian. This will quantitatively justify 
the use of the zero temperature mass for temperatures T < Aqcd- 
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For two light quarks, the chiral Lagrangian, derived from considerations of chiral 
symmetry breaking, is a low energy effective theory of pions, and to lowest order 
(apart from counterterms) is given by 

^2 £2 

C = ^Tr + ^BoTr {MU^ + UM^) , (75) 

where = 96 MeV is the pion decay constant, Bq is a constant describing the scale 
at which chiral symmetry breaking ocurrs, M = diag(m„,mrf) is the two quark mass 
matrix and 

^ z = l,2,3, (76) 

where cxj are the Pauli matrices, tt* are the pions and there is a sum over i. This 
is the lowest order Langrangian and will be all that we will consider here. It is the 
unique Lagrangian satisfying the required symmetries and of dimension two, where 
a derivative has dimension one and the mass matrix has dimension two. There are 
higher order Lagrangians, the next one being of dimension four EH- This Lagrangian 
is only valid for energies where QCD is conhned, and can therefore form bound states 
such as pions. Thus, formulas derived from the chiral Lagrangian are taken to be 
valid only for energies and temperatures below Aqcd- 

To incorporate the axion [32], the mass matrix is modihed by 

M ^ (77) 

where a is the axion and fa is the axion scale. Plugging this and (jZSD into the chiral 
Lagrangian m and expanding one hnds the axion mixing with tts, which, if electric 
charge is incorporated, is the neutral pion. Diagonalizing the mass matrix leads to the 
usual pion mass, plus 0(/7r//a)^ corrections, and the zero temperature QCD axion 
mass, 

ml+ 0{f^/fa)\ (78) 

\JaJ {mu + mdy 

where and md are the up and down quark masses and is the neutral pion 
mass. 

We will calculate temperature corrections using the partition function |33j. In 
particular we will calculate the one loop correction to the scalar potential at hnite 
temperature |33|. The partition function requires moving to imaginary time, inte¬ 
grating imaginary time from 0 to /? = 1/T, where T is the temperature, and taking 
ah helds to be periodic in imaginary time with period /?. 

To evaluate the path integral we make the usual hrst step of shifting the helds by 
their classical solutions. We could expand the chiral Lagrangian (ESI) hrst and then 
make this shift, but it proves easier to parametrize the shift as follows ESI, 

u = a = a + 5, A = AiCJi, i = 1, 2, 3, (79) 

where U and a are the classical helds (i.e. the solutions to the equations of motion) 
and Aj and 6 are the dynamical helds to be integrated over in the partition function. 
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Figure 4: Both curves are plots of the hnite temperature QCD axion mass divided 
by the zero temperature mass. The top line is only valid for T < Aqcd; while 
the bottom line is only valid for T > Aqcd, although each is plotted outside their 
respective domain. 


The partition function leads to an integral that must be done numerically, and is 
part of the one loop corrected mass matrix which must be diagonalized. Upon doing 
so, the axion mass is obtained and can be plotted as a function of temperature, which 
is the top line in figured For the purpose of comparison, also included in the plot as 
the bottom line is the high temperature axion mass, (jH)) or (ESI), which is derived from 
(unconhned) QCD and therefore is only valid for temperatures greater than Aqcd- 
In terms of orders of magnitude, if we take Aqcd rsj O(IO^) MeV, then the top 
line will be valid up to around T ~ 0(10) MeV, where we can see that the axion 
mass barely changes. Where the top line does begin to curve appreciatively, say at 
T ~ 200 MeV, we can no longer trust the plot, as the chiral Lagrangian, from which 
it is derived, is no longer valid. 

Figure m shows that we may take the low temperature axion mass to be equal to 
its zero temperature value for temperatures T < Aqcd- 
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